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ABSTRACT 

Based on a detailed analysis of the high-quality Chandra, XMM-Newton, and 
Suzaku data of the X-ray bright cluster of galaxies Abell 1795, we report clear 
evidence for a two-phase intracluster medium (ICM) structure, which consists of 
a cool (with a temperature T c ps 2.0 — 2.2 keV) and a hot (Tj, ~ 5.0 — 5.7 keV) 
component that coexist and dominate the X-ray emission at least in the central 
80 kpc. A third weak emission component (T 3 ph 0.8 keV) is also detected within 
the innermost 144 kpc and is ascribed to a portion of inter-stellar medium (ISM) 
of the cD galaxy. Deprojected spectral analysis reveals flat radial temperature 
distributions for both the hot phase and cool phase components. These results 
are consistent with the ASCA measurements reported in Xu et al. (1998), and 
resemble the previous findings for the Centaurus cluster (e.g., Takahashi et al. 
2009). By analyzing the emission measure ratio and gas metal abundance maps 
created from the Chandra data, we find that the cool phase component is more 
metal-enriched than the hot phase one in 50 — 100 kpc region, which agrees 
with that found in M87 (Simionescu et al. 2008). The coexistence of the cool 
phase and hot phase ICM cannot be realized by bubble uplifting from active 
galactic nuclei (AGN) alone. Instead, the two-phase ICM properties are better 
reconciled with a cD corona model (Makishima et al. 2001). In this model, the 
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cool phase may be ascribed to the plasmas confined in magnetic loops, which are 
surrounded by the intruding hot phase ICM and have been polluted by metals 
synthesized in the cD galaxy. AGN feedback energy released in the innermost 10 
kpc can serve as the heating source to prevent the loop-interior gas from cooling 
down to temperatures much lower than the observed value. The total gravitating 
mass profile exhibits a hierarchical structure, regardless of the ICM temperature 
modeling. 

Subject headings: galaxies: clusters: individual (Abell 1795) - - intergalactic 
medium — magnetic fields — X-rays: galaxies: clusters 

1. INTRODUCTION 

First revealed by the Einstein observatory, the spectroscopic temperature of the intr- 
acluster medium (ICM) in many cD clusters of galaxies decreases toward the center by a 
factor of 2 — 3 (e.g., Canizares et al. 1988; Mushotzky & Szymkowiak 1988; see also Fabian 
1994 for a review). Generally, the appearance of such a central cool component (hereafter 
CCC) has been described in two possible but competing ways. One is a single-phase (here- 
after IP) scenario, in which a mono-nature plasma showing an inward temperature decrease 
is assumed to permeate the whole cluster (e.g., Allen et al. 2001; Vikhlinin et al. 2006). The 
other is a two-phase (hereafter 2P) scenario, in which the cluster's central region is assumed 
to be occupied by a mixture of two gas components characterized by discrete temperatures 
(i.e., hot and cool), of which the relative volume filling factors slowly vary with radius (e.g., 
Fukazawa et al. 1994, 1998; Takahashi et al. 2009; see Makishima et al. 2001 for a review). 
Although in many cases it is difficult to distinguish between the IP and 2P models due to 
insufficient data quality, the implications of these models on the formation and evolution of 
the CCC are intrinsically different. In the IP scenario the CCC is simply interpreted as the 
dense cluster core, where radiative losses have significantly lowered the gas temperature (see, 
e.g., Peterson & Fabian 2006 for a review), while the 2P scenario favors an interpretation 
that the CCC is caused by the cD galaxy, rather than being a cooling portion of the ICM 
(Makishima et al. 2001 and references therein). Furthermore, the choice between the two 
scenarios has considerable impact on the accuracies of the measurements of both gas metal- 
licities (e.g., the "Fe-bias" discussed in Buote 2000) and, in some extreme cases, gravitating 
mass distributions in clusters. For the latter one, the systematic mass bias induced by dif- 
ferent ICM models can reach ~ 10% for the cluster's central region, which is comparable to 
the typical deviation between the results obtained in X-ray and lensing studies (e.g., Nagai 
et al. 2007; Mahdavi et al. 2008). Hence, it is of fundamental and urgent importance to 
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determine which of the two scenarios is valid, based on the analysis of existing high quality 
X-ray data of nearby, bright clusters. 

So far, only a few works have been done to examine the preference between the IP and 
2P ICM models of cD groups and clusters. Of these, Fukazawa et al. (1994) and Ikebe et 
al. (1999) used the ASCA data to show that a 2P model with temperatures of 1.4 keV and 
4.4 keV for the two components is required to describe the CCC of the Centaurus cluster. 
Their results were recently confirmed by the analysis of the high quality data acquired with 
the European Photon Imaging Camera (EPIC) and Reflection Grating Spectrometer (RGS) 
onboard XMM-Newton (Takahashi et al. 2009). The preference for the 2P scenario was also 
reported in the works on the Fornax cluster (Ikebe et al. 1996; Matsushita et al. 2007) and 
the NGC 5044 group (Buote & Fabian 1998; Buote et al. 2003; Tamura et al. 2003) with 
ASCA, Chandra, XMM-Newton, and Suzaku observations. However, due to the lack of such 
studies on more representative relaxed clusters or groups that exhibit strong CCCs, it is not 
yet clear whether or not the 2P scenario do surpass the IP scenario in general, and if yes, 
to what extent the 2P model can be applied to constrain the origin of the CCC. 

To address these issues, we study in this work the rich, bright cD cluster of galaxies 
Abell 1795 (A1795 hereafter), by jointly analyzing the archival Chandra, XMM-Newton, and 
Suzaku data. A1795 is suited to our purpose as it is a relaxed, nearby (at a redshift of z 
= 0.0625) cluster of galaxies with a luminous X-ray emission (with a 2.0—10.0 keV luminosity 
of Lx ~ 1-0 x 10 45 ergs s -1 ; Xu et al. 1998), which is found peaked at the giant cD galaxy 
PGC 049005 (Jones & Forman 1984; Briel & Henry 1996). Under the 2P assumption, Xu 
et al. (1998) resolved a cool (Tx ~ 1-7 keV) and a hot (Tx ~ 6.5 keV) component in the 
central 220/ifj 1 kpc region of this cluster with ASCA, although Ettori et al. (2002) employed 
the IP model instead to fit Chandra spectra extracted from the same region, and found 
an inward temperature decrease to ~ 2.5 keV. In the past 11 years, X-ray observations of 
this cluster have been accumulated to « 300 ks, 140 ks, and 270 ks with Chandra, XMM- 
Newton, and Suzaku, respectively. Clearly this is an ideal target for utilizing as much of 
existing high-quality X-ray data as possible to judge which spectral model is more correct. 

The layout of this paper is as follows. Section 2 gives a brief description of our data 
reduction procedure. The data analysis and results are described in §3. We discuss the 
physical implication of our analysis in §4, and summarize our work in §5. Throughout the 
paper we assume a Hubble constant of Hq = 71hji km s _1 Mpc" 1 , a flat universe with 
the cosmological parameters of VL M = 0.27 and Q\ = 0.73, and quote errors by the 90% 
confidence level unless stated otherwise. At the redshift of this cluster, 1' corresponds to 
about 72hji kpc. To compare with previous results, we adopt the solar abundance standards 
of Anders & Grevesse (1989). 
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2. OBSERVATION AND DATA REDUCTION 
2.1. X-ray Observation 

iU.i. Chandra 

In Table 1, we list 11 Chandra datasets of A1795 obtained with the advanced CCD 
imaging spectrometer (ACIS), which are all used in the analysis. Among the 11 observations, 
five ACIS-S3 exposures focused on the cluster center, and six ACIS-I exposures covered off- 
center regions up to r ~ 1.5 Mpc. All data were telemetered in VFAINT mode, except 
for the one with ObsID 494 in FAINT mode. Using CIAO v4.4 and CALDB v4.4.7, we 
removed bad pixels and columns, as well as events with ASCA grades 1, 5, and 7. We then 
executed the gain, CTI, and astrometry corrections. The cumulative exposure time was 
reduced from ps 170 ks to ps 163 ks after removing intervals contaminated by occasional 
background flares with count rate >20% of mean value, which were detected by examining 
0.3 — 10.0 keV lightcurves extracted from source free regions near the CCD edges (e.g., Gu et 
al. 2009). When available, ACIS-S1 data were also used to crosscheck the determination of 
the contaminated intervals. The obtained clean exposure time for each observation is listed 
in Table 1. In Figure la we show the combined ACIS image, which has been corrected for 
exposure but not for background. 



2.1.2. XMM-Newton 

XMM-Newton was used to observe the central and the south periphery regions of A1795 
on 2000 June 26 and 2003 January 13, respectively, both with the EPIC operated in the full 
window mode with thin filter, and with the RGS in the spectroscopy mode (Table 1). Data 
reduction and calibration were carried out with SAS vl 1.0.1. In the screening process we 
set FLAG = 0, and kept events with PATTERNS 0-12 for MOS cameras and events with 
PATTERNS 0-4 for pn camera. By examining lightcurves extracted in 10.0 — 14.0 keV 
and 1.0 — 5.0 keV from source free regions, we rejected time intervals affected by hard- and 
soft-band flares, respectively, in which the count rate exceeds a 2a limit above the quiescent 
mean value (e.g., Katayama et al. 2004; Nevalainen et al. 2005). The cleaned MOS1, MOS2, 
and pn datasets of the central pointing have exposure times of 43.6 ks, 40.3 ks, and 30.2 ks, 
respectively (Table 1). The RGS data were screened using the method described in Tamura 
et al. (2001). 
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2.1.3. Suzaku 

As listed in Table 1, A1795 was also observed by Suzaku on 2005 December 10, with 
five pointings aimed at its central region, near south and north regions (offset ~ 12'), and 
far south and north regions (offset ~ 24'). The onboard X-ray Imaging Spectrometer (XIS; 
Koyama et al. 2007) and the Hard X-ray Detector (HXD; Takahashi et al. 2007) were both 
operated in normal modes. The same XIS datasets were already utilized by Bautz et al. 
(2009). In analyzing the obtained data, we used the software HEASoft 6.11.1 and latest 
CALDB (20111109 for the XIS and 20110913 for the HXD). We started with version 2.0 
processing data, and removed the data obtained either near South Atlantic Anomaly or at 
low elevation angles from the Earth rim (< 5° and < 20° for night and day, respectively). 
Cut-off rigidity criteria of > 8 GV for HXD-PIN data were also applied. For the XIS 
instrument, we further examined 0.3 — 10.0 keV lightcurves of a source free region in each 
CCD, and filtered off anomalous time bins with count rate above the 2a limit of the quiescent 
mean value. The combined XIS image is shown in Figure lb. The data obtained in the far 
north field are contaminated seriously by solar wind charge exchange emission (Fujimoto 
et al. 2007), and thus are not included in the following data analysis. The obtained clean 
exposure times are listed in Table 1. We do not consider HXD-GSO data in this paper. 



2.2. Point Sources, Background, and Systematic Uncertainties 

We excluded point sources detected beyond a 3a threshold on the ACIS images with the 
CIAO tool celldetect, and masked the corresponding regions on the EPIC and XIS images 
after considering the differences in the Point Spread Functions (PSFs) of the instruments. 
The mask regions have radii of 12" and 70" for the EPIC and XIS, respectively. The point 
sources outside the ACIS sky coverage were detected and masked using the EPIC images. 
The ancillary response files (ARFs) and the redistribution matrix files (RMFs) were calcu- 
lated with the CIAO tools mkwarf and mkacisrmf for the ACIS, and with the SAS tools 
arf gen and rmf gen for the EPIC. We calculated the RMFs for the RGS using the SAS tool 
rgsrmf gen, after blurring the line spread function by convolving the raw RMFs with the 
surface brightness profile (§3.3.1) extracted from the 0.5 — 8.0 keV ACIS image along the 
dispersion direction. As for the XIS, we generated the response files using xissimarfgen 
(version 2008-04-05) and xisrmfgen (version 2007-05-14). To enhance the precision of the 
Monte-Carlo simulation carried out by xissimarfgen, we created the initial feed photon list 
by combining all the exposure-corrected and background-subtracted ACIS-I images. The 
method to build the response of HXD-PIN is described in §3.4. 

We estimated the background as a combination of three independent components, i.e., 
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non X-ray background (NXB), cosmic X-ray background (CXB), and Galactic emission. For 
all the XIS data, we created the NXB spectra using a dark earth observation database. 
The CXB and Galactic emission components in the XIS data were estimated by analyzing 
a spectrum extracted from a region about 26' - 30' (« 1.9 - 2.2/if/ Mpc at the distance of 
A1795) south of the cluster center, which is covered by the dataset with ObsID = 800012050 
(Table 1); a similar background region was used in Bautz et al. (2009). Since this region 
is outside the virial radius r 2 oo = 1-9/ifj 1 Mpc, where the brightness of cluster emission is 
expected to be ~ 10~ 13 ergs cm -2 s _1 deg -2 in 0.5 — 2.0 keV (Roncarelli et al. 2006), 
lower than the detection limit of the XIS, i.e., ~ 10 -12 ergs cm~ 2 s _1 deg~ 2 (Bautz et al. 
2009), the ICM component can be ignored in the background fitting. After subtracting the 
NXB from the extracted spectrum, we fitted the resulting spectrum with an absorbed power 
law model (photon index T = 1.4) describing the CXB, and two unabsorbed optically thin 
thermal models (abundance = 1 Z & , temperatures = 0.08 keV and 0.2 keV; Snowden et 
al. 1998) describing the Galactic emission. The 2.0 — 10.0 keV CXB flux was estimated as 
6.6 x 10~ 8 ergs cm~ 2 s _1 sr _1 , which is consistent with the results of Kushino et al. (2002). 
The obtained CXB + Galactic background template was applied to all the subsequent XIS 
spectral analysis assuming the same uniform CXB + Galactic emission distribution over the 
field of view (e.g., Sato et al. 2008). 

In order to achieve an accordant background model for the Chandra ACIS, we utilized 
the above CXB + Galactic background template and calculated the ACIS NXB spectra in 
an adaptive way as follows. For all ACIS-I observations, background spectra were extracted 
from regions 17' — 21' (~ 1.3 — 1.5/ifj 1 Mpc) away from the cluster center, where the ICM 
brightness was measured to be weak but cannot be neglected, on a level of 10% to 40% of 
the CXB + Galactic components (Bautz et al. 2009). We subtracted the CXB + Galactic 
components from the extracted spectra, whose brightness was fixed at the XIS result obtained 
above, and fitted the resulting spectra with an empirical NXB model that consists of a broken 
power law and five narrow Gaussian lines (e.g., Humphrey & Buote 2006), plus a thermal 
APEC model for the cluster emission. Since the temperature and metal abundance of this 
ICM component cannot be well constrained, due to the relatively high NXB level of the 
ACIS, we fixed them to previous Suzaku results (temperature = 2.5 keV, abundance = 0.1 
Z Q ) reported in Bautz et al. (2009). In the r = 17' — 21' region, average surface brightness of 
the ICM component was obtained as 1.7 ±0.8 x 10 -12 ergs cm~ 2 s^ 1 deg -2 in 0.5 — 2.0 keV, 
which is about 25% of the CXB + Galactic components, in rough consistent with Bautz et al. 
(2009). The ACIS-I NXB model was thus determined. The same procedure was applied to 
all ACIS-S data, but source free regions on the SI chip were used as the background regions 
instead. To crosscheck above background model with the blank sky templates, we applied 
both backgrounds to the spectral analysis in §3.1.1. By fitting the background-subtracted 
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ACIS spectra extracted from 9'. 8 — 17'. with a single-phase thermal model, we found that 
the ICM temperatures with the two background models differ by 0.35 keV, which is less 
than the la uncertainties of 0.5 keV. The error caused by background model is much smaller 
for the inner regions, where the ICM emission contributes to > 80% of the total counts in 
0.5-10.0 keV (Table 2). 

The background model for the XMM-Newton EPIC data was determined with a method 
similar to that used for the Chandra ACIS data. The EPIC spectra extracted from 19 — 25' 
(~ 1.4 — 1.8h?i Mpc) away from the cluster center, covered by the data with ObsID = 
0109070201, were fit with a model consisting of the CXB + Galactic components, the EPIC 
NXB model including a broken power law and six Gaussian lines (e.g., Gastaldello et al. 
2007; Zhang et al. 2009), and a cluster component approximated by an APEC model. The 
temperature and abundance of the cluster component were fixed at 2.4 keV and 0.1 Z Q 
(Bautz et al. 2009), respectively, and the 0.5 — 2.0 keV surface brightness was obtained 
as 1.4 ± 1.1 x 10~ 12 ergs cm~ 2 s _1 deg~ 2 15% of the CXB + Galactic components), 
in agreement with Bautz et al. (2009) and our Chandra estimate. The best-fit CXB + 
NXB + Galactic models were used to create background spectra. In Figure 2 a, we plot the 
NXB-subtracted background spectra against the best-fit CXB and Galactic components for 
the ACIS, EPIC, and XIS, as well as the cluster emission component for the former two 
instruments. For the RGS data, we employed the blank sky template (e.g., Takahashi et al. 
2009) in the following spectral analysis. 

Errors quoted below in the spectral fittings were estimated by taking into account both 
statistical and systematic uncertainties. The former was calculated by scanning over the 
parameter space with the XSPEC command steppar, as the fitting was repeated for a few 
iterations at each step to ensure that the actual minimum x 2 is found. For the latter, we 
altered the normalization of the CXB spectra by 10% (Kushino et al. 2002) to approximate 
its field-to-field variation, and similarly, re-normalized the NXB components by 2%, 2%, 
and 1% for the ACIS, EPIC, and XIS data (e.g., Hickox & Markevitch 2006; De Luca & 
Molendi 2004; Tawa et al. 2008), respectively, to assess the error ranges. We have also tried 
assigning a systematic error of 2% to approximate the unassessed calibration uncertainties 
of the EPIC data (e.g., Takahashi et al. 2009), which turned out to have negligible impacts 
on the results and thus was not included in the actual fittings. 
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3. ANALYSIS AND RESULTS 
3.1. Azimuthally Averaged Spectral Analysis 

3.1.1. Single-phase ICM Model 

We extracted the ACIS and EPIC spectra from eight thin annuli, i.e., — 30/iyf 1 kpc 
(0'-0'.4), 30-51/if/ kpc (0'.4-0'.7), 51-80/ifj 1 kpc (O'J-l'.l), 80-116/if/ kpc (lM-1',6), 
116-238/^ kpc (l'.6-3'.3), 238-354^ kpc (3'.3-4'.9), 354-707/if/ kpc (4'.9-9'.8) and 
707 — 1335/if/ kpc (9'. 8 — 18'. 5), and the XIS spectra from four thick annuli, i.e., — 144/if/ 
kpc (0'-2'.0), 144-320/if/ kpc (2'.0-4'.4), 320-700^ kpc (4'.4-9'.7) and 700-1444/1^ 
kpc (9'. 7 — 20'. 0). When fitting these spectra, the lower energy cut was set at 0.7 keV, 0.7 
keV, and 0.6 keV for the ACIS, EPIC, and XIS, respectively, while the upper cut was fixed 
at 8.0 keV for all the three detectors. Also, the Si K edge (1.8 — 1.9 keV) was excluded from 
all the XIS spectra. Independent fittings were carried out for the Chandra (ACIS-S and 
ACIS-I), XMM-Newton (MOS and pn), and Suzaku (BI and FI) spectrum sets, by applying 
a common absorbed APEC model and linking all annuli with XSPEC model PROJCT, which 
performs projection of 3-D shells onto 2-D annuli to evaluate the projected emission of outer 
shells on inner ones. For each shell, the gas temperature, metal abundance, and column 
density of the neutral absorber were set free. The best-fit model yielded x 2 / u — 1750/1496, 
1040/881, and 845/775 for the ACIS, EPIC, and XIS spectra, respectively. Figure 2b and 
2c show the best-fit deprojected IP fittings to the spectra extracted from 320 — 700/if/ kpc 
and core regions (0 - 80/if/ kpc for the ACIS and EPIC, and - 144/if/ kpc for the XIS), 
respectively. 

As shown in Figure 3a, the best-fit deprojected IP temperature drops inwards from 
ps 6.0 keV at ps 100 — 350/if/ kpc to ps 3.0 keV in the central 30h^ kpc, which shows 
apparent diagnostic of a cool core (e.g., Sanderson et al. 2006). On the other hand, in the 
cluster's outskirt (~ 350 — 1100/if/ kpc), the temperature declines outwards steeply down 
to ps 3.0 keV. Generally, our IP temperature profiles are consistent with previous reports 
(e.g., Ettori et al. 2002; Vikhlinin et al. 2006; Snowden et al. 2008; Bautz et al. 2009). 
Incidentally, in ps 150 — 300/if/ kpc, the temperature obtained with the ACIS spectra is by 
about 1.0 keV higher than those measured with the EPIC and XIS spectra. This discrepancy 
can be ascribed to a local temperature structure (see §3.2.2), which is smoothed to some 
extent in the XMM-Newton and Suzaku data. 

Although the best-fit IP model can reproduce the three data groups with reasonable 
goodness, i.e., x 2 / u ~ 1-09 ~ 1.18, the fitting residuals become drastically significant in the 
central 80/if/ kpc regions, where the averaged gas temperature is measured to be ~ 4.5 keV 
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(Fig. 3a). As shown in Figure 2 c, the IP model significantly underestimates the Fe-L blend 
and the continuum in > 5.0 keV. This indicates that these annular spectra exhibit stronger 
multi-temperature nature than is predicted by foreground and background contributions 
from different radii that are accounted for by the PROJCT model. To further examine the 
issues with the IP ICM model for these regions, following, e.g., Cavagnolo et al. (2008), 
we fitted the ACIS spectra in 0.7 - 4.0 keV and 4.0 - 8.0 keV with the IP model, both 
extracted from the central 80 /if/ kpc and corrected for projection effects. The IP gas 
temperature obtained in 0.7 — 4.0 keV (T a7 _ 4 . kcV = 3.5 ± 0.2 keV) became lower than 
its counterpart in 4.0 — 8.0 keV (Z4.0-8.0 kcV = 5.7 ± 0.7 keV) at the 99% confidence level. 
The bandpass dependence cannot be ascribed to the calibration problem of Chandra at soft 
band, because a similar dependence can be obtained using the XMM-Newton data, with 
To.7-4.0 kcV — 3.3 ± 0.1 keV and T 4 _8.o kcV — 4.7 ± 0.5 keV. Thus, the IP modeling of the 
ICM is considered inadequate, and a 2P ICM spectral model is hence invoked. 



3.1.2. Two-phase ICM Model 

Next we investigate whether or not the 2P ICM model, as inferred in the above IP 
analysis, is applicable to the spectra extracted in the inner 80/if/ kpc region, after the pro- 
jection effects are taken into account. Here, the 2P ICM components were both represented 
by thermal APEC models, which were constrained to have the same metallicity and ab- 
sorption. Both the cool- and hot-phase ICM temperatures (T c and Th, respectively) were 
tied among all the thin shells in central 80 /if 1 kpc, while the IP ICM model (§3.1.1) was 
retained for the outer parts. This fitting yielded overall x 2 l v = 1608/1494 and 929/879 for 
the ACIS and EPIC spectrum sets, respectively. The fitting goodness was compared to the 
IP one (§3.1.1) using an F-test, which yielded F-statistic of 66.0 and 52.5 for the ACIS and 
EPIC data, respectively, indicating that the 2P model gives a better fit than the IP model 
at > 99% confidence level. As shown in Figure 2c, the 2P ICM model apparently better 
reproduces the Fe-L blend and the continuum in > 5.0 keV than its IP counterpart. This 
shows that the 2P ICM model is strongly and consistently required for the inner 80/if/ kpc 
region even after removing (via PROJCT) foreground and background contributions from 
the outer shells. The 2P temperatures (T c , T h ) were determined as (2.2 ± 0.2, 5.7 ± 0.4) 
keV and (2.0 ± 0.2, 5.0 ± 0.3) keV by using the ACIS and EPIC data, respectively. For the 
regions with r > 80/if 1 kpc, the 2P model does not significantly improve the fitting over the 
IP model. 

The same 2P ICM model was also applied to the thick shell spectra obtained with the 
XIS. Compared to the IP case (§3.1.1; x 2 /" ~ 845/775), the 2P ICM model again gave 
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significantly better fit to the XIS spectra within central 320/if/ kpc region, with reduced 
chi-squared improved to x 2 / u = 802/773. Since the PSF of the XIS is much wider than 
those of the ACIS and EPIC, with the XIS alone we cannot constrain the boundary between 
the IP and 2P regions to a small radius, as we did with the ACIS and EPIC data. To 
examine the PSF effect on the spectral fitting of the innermost thick shell (0 — 144/iyj 1 kpc), 
we performed a ray tracing simulation (e.g., Ishisaki et al. 2007; Reiprich et al. 2009) based 
on the ACIS image, and found that the 0.5 — 8.0 keV counts in this shell are mostly 75%) 
from the emission originated in — 80/if/ kpc, only a small fraction 5%) come from 
116 — 238/ifj 1 kpc region where a hotter component (T x ~ 6.5 keV; Figure 3a) is detected 
with the ACIS data. Hence this hot component cannot bias the 2P result significantly. In 
fact, the best-fit 2P temperatures for this region, T c — 2.1 ± 0.5 keV and Th = 5.5 ± 0.4 keV, 
are consistent with those obtained in the thin shell analysis. Thus, the PSF of the XIS little 
affects the 2P results obtained for the innermost thick shell. 

As the most consistent form of our 2P analysis, we simultaneously fitted the thin shell 
spectra (the ACIS and EPIC) and those from the thick shell (the XIS), by allowing to 
vary independently the pair of 2P temperatures of each inner shell, i.e., < kpc and 

< 320/17/ kpc for the ACIS /EPIC and XIS cases, respectively. Nearly the same fit goodness 
was achieved (Table 2), and the best-fit T c and Th, as shown in Figure 4a and Table 2, exhibit 
nearly insignificant spatial variations across the inner shells. The values of T c and T h are 
consistent among the three instruments, and broadly consistent with the ASCA result, (T c , 
T h ) = (1.7 ± 0.3, 6.5 ± 0.6) keV, reported in Xu et al. (1998). In short, the ICM in the cool 
core of this cluster can be described by two discrete temperatures, T c « 2.0 — 2.2 keV and 
Th ~ 5.0 — 5.7 keV, which are both consistent with being spatially constant. This makes the 
simple 2P ICM picture (Makishima et al. 2001) a natural and reasonable description. 



3.1.3. A Weak 0.8 keV Spectral Component Within the Central 144h 7l kpc 

Employing the IP formalism, Fabian et al. (2001) reported a filamentary structure in 
the inner 50/if/ kpc region of A1795, with a possible temperature of ~ 1 keV, apparently 
lower than the value of T c obtained in our 2P analysis. To look for such a component, we 
added a third APEC component to the 2P ICM model describing the XIS data from the 
central 144/iyj 1 kpc region. Indeed, we obtained a sig nificantly better fit (x 2 /v = 791/771) 
by adding a third APEC component, whose temperature is kT 3 = 0.8±0.4 keV, while its 
metal abundance, absorption, and redshift were tied to those of the 2P ICM components. 
An F-test indicates that the probability of this improvement being caused by chance is 
< 5 x 10~ 3 . The 0.3 - 10.0 keV luminosity of this 0.8 keV component is 3.6+2:6 x 1q42 er S s 
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s 1 , which is consistent with the luminosity of the filamentary structure measured with the 
ACIS in Fabian et al. (2001; ~ 4 x 10 42 ergs s" 1 ). As shown in Table 2, the values of T c 
and Th, as well as the ICM abundances, remain nearly the same by adding this 0.8 keV 
component, because it contributes only m 0.7% of the total 0.3 — 10.0 keV luminosity of the 
central 144/if/ kpc region 5.0 x 10 44 ergs s _1 ). 

In order to crosscheck the XIS result, we analyzed the XMM-Newton RGS and depro- 
jected EPIC spectra extracted from the central 80/if/ kpc region by fitting them with three 
component spectral model, i.e., 2P ICM plus a third APEC component. Since the RGS is 
not so sensitive to hot gas components, we fixed T c and Th at the best-fit EPIC results, i.e., 
2.0 keV and 5.0 keV, respectively. By fitting the RGS1 and RGS2 spectra simultaneously 
in 6 — 23 A, a weak 0.8 keV component was again detected, with a probability of 6 x 10~ 2 
for the detection to be caused by chance. The fitting of the EPIC spectra was improved 
to x 2 /V = 923/878 by including the 0.8 keV component, significantly better than previous 
2P fitting in terms of F-test (> 98% confidence level). The 0.3 — 10.0 keV luminosities of 
this component determined with the RGS and the EPIC are 2.7 ± 2.3 x 10 42 ergs s" 1 and 
2.6 ± 1.7 x 10 42 ergs s _1 , respectively, nicely consistent with the XIS value and the ACIS 
result in Fabian et al. (2001). 

3.1. 4- Filling Factor of the Cool Phase 

Using the high resolution Chandra, XMM-Newton, and Suzaku data, we have demon- 
strated that the 2P ICM model gives a significantly better description of the ICM thermal 
condition in the central 80/if/ kpc of A1795 than the IP counterpart. The cool and hot 
phase ICM components have 0.3 — 10.0 keV luminosities of about 1.4 x 10 44 ergs s -1 and 
2.5 x 10 44 ergs s _1 , respectively. Also a weak 0.8 keV component has been detected in the 
same region, which has a 0.3 — 10.0 keV luminosity of about 3.6 x 10 42 ergs s -1 . 

The 2P ICM picture implicitly assumes that these phases with different temperatures 
coexist, each occupying a certain fraction of the total volume under study in the cluster 
center, which can be described as volume filling factor (e.g., Ikebe et al. 1999; Makishima 
et al. 2001; Takahashi et al. 2009). To quantify this quantity, we followed the method 
described in Ikebe et al. (1999) and calculated the filling factor of the cool phase gas, r] c (R), 
where R denotes 3-D radius. With i] c , the specific emission measures of the two phases, 
Q C (R) and Qh(R), are described as 

Q C (R) = n c (R)\(R), Q h (R) = n h (R) 2 (l - Vc (R)), (1) 

where n c (R) and rih(-R) are the density distributions of the cool and hot phases, respectively. 
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Assuming a pressure balance between the two phases 



n c (R)T c (R) = n h (R)T h (R), 



(2) 



we have 



Vc(R) 



l + 



T h (R)\ ( Q h (R) 



T C (R)J \Q C (R) 



(3) 



Figure 46 shows radial profiles of rj c (R), calculated with the best-fit deprojected ACIS, 
EPIC, and XIS model parameters as listed in Table 2. Thus, the ACIS and EPIC datasets 
consistently indicate that the hot component dominates in volume over its cool counterpart 
not only in outer regions, but also in the core region (< kpc). The cool phase gas 

occupies up to r/ c ps 0.2 of the volume in the central 30/ifj 1 kpc, whereas r/ c declines steeply 
to below 0.05 outside the central 80/if/ kpc region. The XIS result for the inner 144/iyf 1 
kpc is consistent with those of the ACIS and EPIC, whereas the relatively large value of 
rj c indicated by the outer 144 — 320/if/ kpc XIS bin can be attributed to the broad PSF of 
the XIS. In fact, by performing a ray tracing simulation (see §3.1.2 for more details), the 
outer XIS bin at 144 - 320/if/ kpc was found to be contaminated by photons scattered from 
the inner regions, in such a way that 37% of the emission in this region actually comes 
from the central 144/iy^ kpc. In this case, a weighted mean of r/ c ps 0.08 measured with 
the XIS over inner 144/z,^ 1 kpc and r/ c ps 0.01 measured with the other two missions at the 
144 - 320/ifi 1 kpc region becomes ~ 0.04, in agreement with the outer XIS point. Hence, 
after correcting for the PSF effect, all data indicate the same 2P configuration in the central 
region. 



3.1.5. 2P vs. Multi-phase ICM Model 

The ICM in the central 144/ifj 1 kpc may alternatively in a multi-phase condition (e.g., 
Kaastra et al. 2004), to be described by a more complicated emission measure distribution 
in a wide temperature range than that used in §3.1.2. To examine this possibility, we 
analyzed the deprojected XIS spectra from the inner 144/ii^ 1 kpc region with the same multi- 
temperature fitting approach as used in Tamura et al. (2001) and Takahashi et al. (2009). 
That is, the spectra were modeled as a cumulative contribution of seven APEC components, 
whose temperatures are given as T , 1.5T , (1.5) 2 T ,..., and (1.5) 6 T , where T is a base 
temperature and left free in the fitting. The seven components were constrained to have 
the same metal abundance, and suffer the same absorption. The fit has been acceptable, 
and the goodness of fitting (x 2 / u = 786/769) is as good as that of the 2P ICM plus a 0.8 
keV component model (Table 2). The base temperature was constrained as T = 0.7 ± 0.2 
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keV. As shown in Figure 5, in the multi-phase model, only one weak component at ~ 0.7 
keV and two strong components at 2.4 keV and 5.3 keV remain significant (68% confidence 
level), while the rest components cannot be constrained. This is essentially identical to the 
2P ICM plus 0.8 keV modeling. Our result is consistent with the multi-phase fitting with 
the EPIC data in Kaastra et al. (2004), which also shows a two-temperature structure (2.5 
keV and 4.9 keV, see their Table 6) in the cluster center. Hence, the XIS data prefers the 
discrete 2P model, to a continuous temperature distribution, although the latter cannot be 
ruled out based on available data. 



3.1.6. Metal Abundance and Absorption Distributions 

As shown in Figure 3b, the deprojected metal abundance profiles appear roughly consis- 
tent between the IP and 2P ICM modelings (§3.1.1 and §3.1.2, respectively). Both profiles 
show a peak in the shell of 30 — 51/if/ kpc, and a mild decline outwards. This abundance 
profile, however, should be regarded as an average among those of different elements, because 
we have so far assumed the solar abundance ratios. To examine the ICM for possible devi- 
ation from the solar ratios, we employed two VAPEC models and reran the deprojected 2P 
fittings of the Suzaku XIS spectra extracted from the central 320/if/ kpc region. Specifically, 
we left the abundances of O, Mg, Si and Fe free, fixed the abundances of He, C and N at the 
solar value (Anders & Grevesse 1989), and tied the abundance of Ni to that of Fe and those 
of other elements (Ne, Al, S, Ar and Ca) to that of Si. This model gave nearly the same set 
of temperatures for the 2P ICM, and a fit goodness (x 2 /^ = 787/761) slightly better than 
that of the best-fit APEC model {^jv = 802/773). As shown in Table 3, the best-fit Fe and 
Si abundances increase significantly towards the center, while the O and Mg abundances are 
nearly constant throughout the cluster. This confirms the previous ACIS result reported in 
Ettori et al. (2002; see also Matsushita et al. 2007). We also added another APEC model (A 
= 0.5 Z Q ; see Table 2) in the central 144/i^j 1 kpc to account for the 0.8 keV component, while 
the model gave nearly the same best-fit abundance profiles for the 2P ICM components. 

To examine the possible existence of any intrinsic absorption, we compared the ab- 
sorption obtained in the spectral analysis with the Galactic value. Regardless of the IP or 
2P modeling, the ACIS and XIS fitting results revealed a significant excess absorption by 
1 —2 x 10 20 cm -2 (Table 3) beyond the Galactic value of 1.2 x 10 20 cm" 2 , at least in the central 
30/17/ kpc region. The excess absorption could not be mitigated by varying the abundances 
of specific elements (e.g., oxygen) whose emission lines couple with the absorption feature. 
This result is in good agreement with the ACIS result of Ettori et al. (2002; ps 2.5 x 10 20 
cm -2 ). On the contrary, our best-fit EPIC result 1.1 x 10 20 cm -2 ), as shown in Table 3, 
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implies lack of excess absorption, in agreement with Nevalainen et al. (2007) who used the 
same XMM-Newton MOS data. The discordance among the three detectors maybe caused 
by a temporal solar wind charge exchange emission, which biased the EPIC absorption low. 
As shown in Wargelin et al. (2004), the brightness of charge exchange can reach ~ 2 x 10~ 6 
photon s _1 arcmin -2 cm -2 in 0.5 — 0.9 keV, which is sufficient to undermine the absorption 
by 0.5 — 1.0 x 10 20 cm~ 2 . The origin of the central excess absorption detected with the ACIS 
and XIS, on the other hand, still remains unclear. 

3.2. Projected 2-D Spectral Analysis 

As shown above, the azimuthally-averaged spectral analysis prefers a view of the 2P ICM 
in the central 80/if/ kpc of A 1795 than the IP counterpart. However, we cannot exclude at 
present the possibility that the 2P preference is artificially caused by an anisotropic spatial 
distribution of IP gas temperature that fluctuates between the range from T c to T h . To 
address this issue, we examine two-dimensional (2-D) temperature distribution in the central 
region of A1795. Below, 2-D ICM abundance and Q c /Qh distributions are also presented. 

3.2.1. Analysis Procedure 

Combining five Chandra ACIS-S pointings onto the central region together, we have col- 
lected sufficient photons for a high resolution 2-D spectral analysis. Following the procedure 
described in detail in Gu et al. (2009), a set of > 10000 discrete points (r ; , i = 1, 2, 3, ...), or 
"knots" , were chosen in the central 240/if/ kpc, which are randomly distributed with a sep- 
aration of < 6h^i kpc between any two adjacent knots. To each knot we assigned a circular 
region, which is centered on the knot and has an adaptive radius of 15 — 80/ifj 1 kpc, ensuring 
that it encloses > 10000 photons in 0.5 — 8.0 keV after all detected point sources (§2.2) 
were excluded. The spectrum extracted from each circular region was fitted with the IP 
and 2P APEC models, both subjected to an absorption that was set free. In the 2P spectral 
fitting, the temperatures of the cool and hot components were fixed at 2.2 keV and 5.7 keV 
(i.e., best-fit ACIS results; §3.1.2), respectively, and the two components were assumed to 
have the same abundance and absorption. For each knot r ; , we obtained the best-fit IP gas 
temperature, 1P/2P abundance, and specific emission measure ratio between the cool and 
hot components (Q c /Qh), as well as their la errors. 

Then, following Gu et al. (2009), we calculated continuous maps based on the obtained 
knots. For any position r within the map region, we defined a scale s(r), so that there 
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are > 10000 net photons in a circular region centered at r, whose radius is s(r). The IP 
temperature at r was then calculated by a weighted mean of all the knots in the circular 
region, 

T(r) = ^(Gr^Mv^/J^Gr.ARrs.X when R r , ri < s(r), (4) 

where R rjr . is the distance from r to r ; , and G r . is the Gaussian kernel whose scale parameter 
a is fixed at s{v\). The use of compact Gaussian kernel guarantees an angular resolution 
of ~ kpc within central 100/if/ kpc region. The abundance and Q c /Qh maps, along 

with their la error maps, are calculated in the same way. The resulting IP temperature and 
abundance maps are shown in Figure 6, and the 2P abundance (hereafter A^p) and Q c /Qh 
maps are presented in Figure 7. The 0.8 keV component (§3.1.3) was not considered in the 
above 2-D analysis, because it will introduce negligible effects, i.e., uncertainties of about 
0.05 keV and 0.02 Z & , to the temperature and abundance measurements, respectively. 



3.2.2. Relaxed Cool Core and SE High Temperature Arc 

The primary purpose of 2-D spectral analysis is to assess the validity of spherically sym- 
metric temperature distribution that we assumed in the radial 1P/2P analysis. To do this, 
we examined the obtained IP temperature map for any significant anisotropic distribution. 
As shown in Figure 6 a, the IP temperature distribution in the cluster central region exhibits 
an approximate elliptical symmetry; the temperature variation in the azimuthal direction is 
about 0.5 keV, 1.0 keV, and 1.0 keV at r = 20/if/ kpc, 40/if/ kpc, and 80/if/ kpc, respec- 
tively, which is insufficient, at least by a factor of three, to imitate the obtained 2P ICM 
condition. Similar morphology is seen in the Q c /Qh map of the core region as shown in 
Figure 7b. This confirms that the 2P view is not an artifact caused by a IP ICM with large 
azimuthal asymmetry. It also indicates that the two phases, in terms of the 2P view, must 
be separated on scales smaller than the spatial resolution allowed by the present analysis 
(?s 15/if/ kpc). 

The most prominent feature on the IP temperature map (Fig. 6 a) is a high-temperature 
arc 7.5 keV according to the IP model) located at ~ 100 — 180/if/ kpc southeast (SE) of 
the cluster's center, with an open angle of >120°. Given the temperature difference of about 
2 keV, the hot structure is significant over the ambient on 90% confidence level. This feature 
agrees with a high temperature bin (116 — 238/if/ kpc) on the deprojected ACIS temperature 
profile shown in Figure 3 a. A same high temperature structure was found by Markevitch 
ct al. (2001), who carried out a IP azimuthal spectral analysis for the southern part of 
the cluster. Markevitch et al. (2001) also reported a density jump by a factor of 1.3 — 1.5 
near the inner edge (~ 85/if/ kpc from the cluster center) of the high-temperature arc, and 
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ascribed the density jump to a cold front. Since both the high-temperature component and 
the cold front locate outsides of the 2P region (i.e., r < 80h^ kpc), they are unable to 
affect the obtained 2P result significantly, even for the XIS thick shell as indicated by the 
ray tracing simulation shown in §3.1.2. 



3.2.3. A Possible Correlation Between Metal-rich and Cool-phase Gas 

A comparison of the A 2 p map (Fig. 7a) with the Q c /Qh map (Fig. 7b) suggests a spatial 
correlation on a scale of 50 — 100/if/ kpc between the two quantities. Both maps reveal 
strongly inhomogeneous distributions outside the relaxed core region, with substructures 
protruding towards southwest, northeast, and northwest of the cluster center by ~ 80 — 
lOO/if-, 1 kpc. Most of these substructures show both higher metal abundances, typically by 
a factor of 2 — 3, and larger cool phase fractions, than their neighborhoods. The best- 
fit A 2 p and Q c /Qh values for all the knot-centered circular regions in 50 - 100/if/ kpc, 
obtained in our 2-D spectral analysis (§3.2.1), were compared directly in Figure 7 c. Indeed, 
it reveals a positive linear correlation, which is represented by an analytic form as A 2 p = 
0- 38 -(U2 x Qc/Qh + 0.48±o;q7 z &- The linear correlation coefficient was obtained as 0.97. 

To quantify the suggested correlation, following, e.g., Shibata et al. (2001), we calculated 
a 2-D cross correlation function as, 

/ {A 2P ( ri ) - A^} {g c (r 2 )/Q h (r 2 ) - O/ft) 

= ( -L— L 

\ A 2F Q C /Q h 

where A 2 p and Q c /Qh are averages of A 2 p and Q c /Qh maps, respectively, r\ 2 = \r\ — r 2 \ is the 
distance between ri and r 2 , and the bracket is ensemble average. The obtained correlation 
function is shown in Figure 7d. A reference profile (£ re f(R)) was also calculated, based on a 
series of random maps, A 2 p iran and Q c ,ran/Qh,ran, which were obtained by randomizing the 
maps within the same data range and smoothing to the same spatial resolution as the original 
A 2 p and Q c /Qh maps. The error bars shown in Figure 7<iwere calculated from the variances 
of £(R) when scattering the A 2P and Q c /Qh maps by their la error maps. Comparing to 
the £ re f(R), the £(R) profile shows a significant excess within lOO/i^ 1 kpc, reconfirming the 
result that cool-phase ICM is indeed more metal-rich. 

Given the detected correlation, the assumption made in §3.1.2 and §3.1.6 that the two 
phases ICM have approximately the same abundances may not stand now. Therefore, we 
refitted the deprojected XIS spectra of the central 144/iy/ kpc region with the 2P ICM 
(VAPEC + VAPEC) + 0.8 keV component (APEC; 0.5 Z G ) model, same as the one in 
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§3.1.6, except that the iron abundances of the cool and hot phases were let float separately. 
The absorption column densities of the three components were again tied together in the 
fitting as a free parameter. We found that the abundance of cool phase ICM (Ap efi ) does 
appear higher than that of the hot one (Ap e ,h); i-e., (^4Fe,c, ^Fe,h) — (0.80 ± 0.25, 0.36 ± 0.06) 
Zq. Compared to the model presented in §3.1.6 (Table 3), the current model improves fit 
goodness to \ 2 1 v — 779/759 from 787/761, with a probability of 2 x 10~ 2 for the improvement 
to be caused by chance. 

Following, e.g., Simionescu et al. (2008), the iron abundance to be obtained with a 2P 
model assuming a single common metallicity can be estimated by 

£ ^ Qc^Fc,c + Qh^Fc,h ^ 

Eliminating Q c and Qh with Eq.(l), we have 

a> „ T* Vc A Fe , c + T?(l- Vc )A Fc , h 

Adopting f] c = 0.08 (Fig. 4b) and (Ap SjC , Ap e ,h) obtained above, A' Fe is then calculated as 
0.51Z©, which agrees with that derived with the XIS spectra (0.48±0.03 Zq, Table 3). Thus, 
all the ACIS and XIS results obtained so far can be interpreted as evidence for the relatively 
high metal abundance of the cool phase ICM. 



3.3. Hierarchical Gravitational Potential Structure 

3.3.1. Central Excess in X-ray Surface Brightness 

In many cool-core clusters, the X-ray surface brightness exhibits an central excess over a 
/3 model that well fits the outer region. Makishima et al. (2001) argued that the central excess 
is a combined result of two major effects; the presence of a cool phase component, and the 
existence of hierarchical potential structure. According to their definition, the hierarchical 
potential is specified as a halo-in-halo structure, i.e., a smaller potential component is nested 
on a larger one, so that the overall potential shows a central deepening relative to a simple 
King profile. Here we examined the surface brightness profile of A1795 for such a central 
excess. In order to measure the surface brightness profile precisely by taking advantages of 
both high spatial resolution of the Chandra ACIS and stable low background of the Suzaku 
XIS, we calculated exposure-corrected surface brightness profiles from the inner 1000/if/ kpc 
ACIS and 1200 — 1800/if/ kpc XIS data, i.e., Sacis(t") an d •S'xis(^), respectively, where r is 
projected radius. The method described in, e.g., Markevitch et al. (1998), was employed 
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to compensate discrepancies on calibration and instrumental background between the two 
instruments. First, we performed simulations using the xissim tool (Ishisaki et al. 2007) 
to smooth the Chandra images with the Suzaku PSF, and extracted the simulated surface 
brightness profile in the central lOOO/i^ 1 kpc, i.e., S' Acls (r). Then, S' ACIS (r) was scaled 
to Sxis( r ) by solving aS' Acls (r) — b = Sxis( r )> where a and b represent the differences in 
normalization and NXB, respectively. A modified S'acis(^) with the original ACIS resolution 
was obtained by applying a and b to SacisM m the same way as we normalized S' ACls (r). 

Then, the modified S , Acis( r ) an d 5xis( r )> shown in Figure 8, were fitted jointly with 
the (3 model, S(r) = S(0){1 + (r/r c ) 2 }~ 3/3+1 / 2 + Sbg, where 5*(0) is the central brightness, 
r c is the core radius, and Sbg is the background. To represent the possible halo-in-halo 
structure, a double-/? model was also tested as S(r) = ,Si(0){H- (r /r c i) 2 }~ 3/3l+1/2 + S 2 (0){1 + 
(r/r c2 ) 2 }~ 3/32+1 / 2 + S , BG, where subscripts 1 and 2 denote compact and extended components, 
respectively. As shown in Table 4, the fit with the f3 model can be rejected on the 95% 
confidence level, because it significantly underestimates the surface brightness in central 
~ 100/if/ kpc, while the double-/? model gives an acceptable fit to the data. The successful 
fits in Figure 8 are based on a still more sophisticated model, a /3+double-/3 gas density 
model, to be explained later. The properties of the central excess brightness, represented by 
the compact (3 component ((3\ ~ 0.75, r cl 53/if/ kpc), are consistent within errors between 
the 0.5 — 3.0 keV and 3.0 — 8.0 keV bands. Such an energy-independent central excess is 
likely to be caused by an intrinsic hierarchical potential shape, rather than the presence of 
cool phase component. This result agrees with the earlier ASCA result reported in Xu et al. 
(1998). 

To further quantify the central emission excess, we corrected the surface brightness pro- 
file for the projection effect using the best-fit deprojected 2P gas temperature and abundance 
profiles (Table 2; §3.1.2). The observed surface brightness profile was modeled as 



S(r) 



A(T C ,A C )Q C (R) 



RdR 



+ 



A(T h ,A h )Q h (R) 



RdR 



+ 5bg, (8) 



where A is the cooling function, Q c and Qh are the same as in Eq.(l), and Sbc is the averaged 
background value. Following Ikebe et al. (1999), we represented the specific emission measure 
profiles by a 2P /3+double-/3 model, which consists of a ft component for the cool phase ICM, 



Qc(R) = 




-3/9 c 



R < 80hj{ kpc, 
R > 80/if/ kpc, 



(9) 



and a double-/? component for the hot phase ICM, 



Qh(R) — <5o,hi 



R 



R 



c,hl 



-3/?hi 



Qo 



h2 



R 



R 



c,h2 



-3/3 h 



(10) 
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where Qo, c , Qo,hi, and Qo,h2 are the model normalizations of the cool phase, the compact 
hot phase, and the extended hot phase components, respectively. The filling factor r] c is 
renormalized into these parameters. 

With these preparations, we fitted the 0.5 — 8.0 keV ACIS + XIS surface brightness 
profiles using Eq.(8), where A was calculated from the best-fit 2P spectral parameters. The 
fit has been successful with x 2 / u — 161/151. The best-fit values of R cc , R c ,hi, R c ,h2, Pa 
Phi, and /?h2 are given in Table 4. The fittings with surface brightness profiles extracted in 
different energies (i.e., 0.5 — 3.0 keV and 3.0 — 8.0 keV bands) are shown in Figure 8. As 
presented in Table 4, all the fittings yielded consistent parameters for the 2P /3+double-/S 
model. 

In some clusters, the central excess can be alternatively explained by a cuspy dark 
matter distribution proposed in Navarro, Frenk, & White (1996; hereafter NFW model). 
The NFW model predicts potential structure with a single spatial scale, instead of the dual 
structure assumed by the /3+double-/9 model. Given the NFW dark matter density profile, 

Pbm(R) = PcS c + , (11) 

where p c is the critical density, 5 C is characteristic density, and R s is scale radius, the ICM 
density profile is expressed as 

+ R - 1 

n NFW (i?) — n ,NFW it ) (12) 

e«s — 1 

where n ,NFW is model normalization and a is a parameter related to the ICM temperature. 
Then the projected ICM surface brightness profile was modeled as a single component, 

/oo RdR 
A(T 1P ,A lP )n 2 ^ FW (R)-^=== + S BG . (13) 

We found this NFW model cannot give acceptable fits to the 0.5 — 8.0 keV surface brightness 
profile for the entire cluster, with a minimum x 2 / v — 307/155. As shown in Figure 8, the 
central excess emission is too strong to be explained by the NFW model that best-fits the 
> 100/ifj 1 kpc region. Hence, we reconfirmed the ASCA result reported in Xu et al. (1998) 
that the halo-in-halo hierarchical model is preferred to the single-halo NFW model in A1795. 

3.3.2. Central Excess in Total Gravitating Mass 

Next we calculated the total gravitating mass profile M(R) based on the best-fit 2P tem- 
peratures (§3.1.2) and the emission measure profiles obtained with the /3+double-/3 model 
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(§3.3.1). With the help of Eqs.(l)-(3), the best-fit emission measure profiles, Q C (R) and 
Qh(R), were converted to gas density profiles of the cool and hot phases, n c (R) and nh(R), 
respectively. Then, under the assumptions of spherical symmetry and hydrostatic equilib- 
rium, the gravitating mass profile M(R) was calculated as 



M(R) = 



R 2 dP(R) 
Gp g (R) dR ' 



(14) 



where P(R) = n c (R)kBT c (R) = n h (R)kp,Ti l (R) is the gas pressure, p g (R) = pm p {r] c (R)n c (R)+ 
[1 — r) c (R)]nh(R)} is the averaged gas mass density, p = 0.609 is the approximated mean 
molecular weight, and m p is the proton mass. As shown in Figure 9, the obtained mass 
profile exhibits a significant shoulder-like feature at R x ~ 120/if^ kpc, enclosing an excess 
mass (M cxccss ) of about 1.5 x 1O 13 M above a flat core mass distribution calculated with the 
extended hot phase component alone. Our result agrees with the ASCA result, R x « 110/if/ 
kpc and M exce ss ~ 2.1 x 10 13 M o , as reported in Xu et al. (1998). 

To examine systematic errors of the mass profile due to different modelings of the ICM 
temperature distribution, we also calculated the gravitating mass profile using the best-fit 
deprojected IP spectral parameters (§3.1.1). The projected X-ray surface brightness profile 
was modeled with a single ICM component as 



A(T 1P , A 1P )Q 1P (R)-^^= + S BG , 



(15) 



where the specific emission measure is given by a double-/? model, 

-3/3i 
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(16) 



The fit is as good as the 2P case, with \ 2 l v — 165/154. The resulting mass profile, derived by 
substituting P(R) = nip(R)kp,Tip(R) and p g (R) = pm p nip(R) to Eq.(14), is shown in Figure 
9 by a red curve. Thus, the shoulder-like potential structure is again found, although the 
IP model systematically underestimates the gravitating mass within 80/if/ kpc by rs 30%. 
This bias is negligible in > 80 kpc. 

Another concern is that, the assumption of hydrostatic equilibrium may not hold for 
the ICM near the cold front located southeast of the cluster center (§3.2.2), casting doubt 
on the validity of mass profile obtained in the related region (Markevitch et al. 2001). To 
investigate this, we excluded the southern half of the cluster and recalculated mass profile 
for the northern half using the best-fit 2P spectral results and surface brightness profile 
with the ACIS data. The 2P temperatures, T Cinort h = 2.0 keV and 7h,north = 5.9 keV, are 
consistent with previous ACIS results. The same is the obtained nh,north(-R) profile, except 
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for the annulus of the cold front (i.e., R = 100 — 180/i 71 kpc), where rih,north(-R) is lower than 
n h {R) by 5 - 20% in 100 - 180/iff 1 kpc. After comparing M nort t 1 (-R) to M(R) in Figure 9, we 
found that the shoulder-like structure, though diminished by w 10%, still remains significant 
after the high-temperature arc is excluded. Hence we conclude that the existence of central 
hierarchical structure is unambiguously confirmed in A1795. 

3.4. Hard X-ray Emission Component 

Here we examined the 12.0 — 50.0 keV HXD-PIN data for possible existence of any 
extra emission component. First, following, e.g., Nakazawa et al. (2009), we calculated the 
effective area of the HXD-PIN considering the location and extension of the source. Based 
on the observed best-fit gas temperature and gas density profiles (Fig. 4 a and Table 4), 
we calculated the ICM surface brightness distribution map Shard^, y) in 12.0 — 50.0 keV. 
Then we divided the whole HXD field of view into 1' x V blocks, calculated the point source 
ARF for each block at (x, y) using the ftool hxdarfgen, and convolved all the monochromatic 
ARFs with the normalized Shard^, y)- The PIN/XIS cross normalization factor of 1.132 
(Suzaku Memo 2007-110) was also adopted in the resulting ARF. To assess the NXB level, 
we utilized the "tuned" NXB model, which provides an optimized reproductivity by making 
use of HXD-GSO information (Suzaku Memo 2008-03@). The CXB model was assumed as 
N(E) = 8.69 x 10" 4 x (E/1.0 keV)" L29 x exp(-£/40.0 keV) photons cm" 2 s" 1 keV -1 , where 
the normalization was set to match the HXD-PIN opening angle of 0.32 deg 2 (e.g., Nakazawa 
et al. 2009). 

As shown in Figure 10, the NXB-subtracted XIS and HXD-PIN spectra in 4 — 50 keV 
were tentatively fitted with a CXB + ICM model. The CXB component was modeled by a 
cutoff powerlaw with index of 1.29 (Boldt 1987), and the ICM component was represented 
by the 2P model derived in §3.1.2. The 2P temperatures were set to (T c , Th) = (2.1, 5.5) 
keV, and their relative normalizations were fixed at the best-fit value in our 2P analysis 
with the XIS. The fitting is acceptable with x 2 / u = 45/55, and no additional hard X-ray 
component is required in 12 — 50 keV band. To obtain an upper limit on any hard X-ray 
component, we fitted the combined XIS-HXD spectra with a 2P ICM plus T = 2 power-law 
model (Nakazawa et al. 2009). Taking account of the NXB systematic uncertainty of 2.0% 
at 90% confidence level (Fukazawa et al. 2009), the upper limit flux of power-law emission 
in 12.0 — 50.0 keV band was estimated to be 8.2 x 10~ 12 ergs cm -2 s _1 . A similar upper limit 



1 http://www.astro.isas. ac.jp/suzaku/doc/suzakumemo/suzakumemo-2007-ll.pdf 
2 http://www. astro, isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo-2008-03.pdf 
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on any thermal hard X-ray excess, 6.2 x 10 12 ergs cm 2 s 1 , was obtained by alternatively 
fitting the combined spectra with 2P ICM plus a 10 keV APEC model. 

4. DISCUSSION 

We have analyzed Chandra, XMM-Newton, and Suzaku data of the relaxed galaxy clus- 
ter A 1795. The deprojected spectra extracted from the central 80 h^l kpc region have been 
successfully reproduced by invoking two major ICM components with discrete temperatures, 
indicating a 2P property. A third weak 0.8 keV component is marginally detected in the core 
region. In §3.2, we analyzed the Chandra data, under 2-D (spatially) and 2P (spectroscop- 
ically) formalism. The ICM metallicity was found to show a similar spatial distribution to 
that of the cool phase component, indicating that the cool phase ICM is more metal-enriched 
than the hot phase one. 

4.1. Two-phase ICM Properties 

A joint fitting of the thin shell spectra from the ACIS and EPIC, and thick shell spectra 
from the XIS, consistently indicates a clear preference of the two-phase ICM nature, over 
the single-phase one, in the central 80 kpc region of A1795. The ICM therein can be 
characterized by two representative temperatures, T h = 5.0 — 5.7 keV and T c = 2.0 — 2.2 keV 
for the hot and cool phases, respectively, whose 0.3 — 10.0 keV luminosities were estimated to 
be L h = 2.5 x 10 44 ergs s -1 and L c = 1.4 x 10 44 ergs s _1 . Both components show insignificant 
temperature gradients in the 80 kpc region. The hot phase ICM dominates in volume 
over the cool phase one by a factor of > 4, even in the innermost region. The 2P view is 
consistent with previous ASCA results reported in Xu et al. (1998). Employing the 2P ICM 
model and fitting the cool phase and hot phase components with /3 and double-/? gas density 
models, respectively, we quantified the shape of cluster potential. As shown in Figure 9, 
the total gravitating mass profile exhibits a central mass excess (M cxccss ^ 1.5 x 1O 13 M ), 
which has a similar spatial scale to the region of clear 2P property. This indicates a potential 
hierarchy, with a lOO/i^ 1 kpc level halo nested in the center of a cluster-scale one. 

As shown in Figure 3b, the ICM metallicity is also enhanced in the 2P region. Fur- 
thermore, by comparing the 2-D metallicity and the emission measure ratio maps created 
with the 2P ICM model, we found a possible spatial correlation on > la level between the 
cool-phase fraction (or rj c ) and the ICM metallicity in the 50 — 100/if/ kpc region. The 
2-D maps also revealed a filamentary distribution of the cool, metal-rich gas in this region. 
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The subsequent 2P fitting to the XIS spectra confirmed the metallicity enhancement in the 
cool phase relative to the hot phase, giving the best-fit iron abundances of the cool phase 
and hot phase ICM to be 0.80 Z & and 0.36 Z Q , respectively. Such a picture resembles the 
previous findings of the metal-rich, multi-phase arms in M87 (Simionescu et al. 2008). In 
the innermost core of A1795, i.e., r < 25/if/ kpc, however, a dip is seen in the abundance 
distribution, which does not correlate with the cool phase fraction. Such a feature is also 
seen in the coolest spots of some other clusters, while its origin still remains unclear (e.g., 
Sanders et al. 2006). 

In §3.1.3 and §3.1.5, we showed that the XIS and RGS spectra of the central 144/iyf 1 
kpc region can be better reproduced by including an additional weak 0.8 keV component, 
without significantly revising the properties of the cool and hot ICM components. This 0.8 
keV component has a 0.3 — 10.0 keV luminosity of 3.6+2% x 10 42 er S s s_1 > an d is dimmer than 
the hot ICM component by two orders of magnitude. Since the temperature and luminosity 
of this component are consistent with those of an X-ray/Ho; filament in the core of A1795 
as reported in Fabian et al. (2001), the 0.8 keV emission might be associated with a portion 
of ISM of the cD galaxy, which is presumably confined within some filamentary magnetic 
structures. 

The 2P formalism has also been successful for the central region of the Centaurus cluster, 
as firstly reported by Fukazawa et al. (1994) and Ikebe et al. (1999) with ASCA data, and 
recently reinforced by Takahashi et al. (2009) using XMM-Newton. According to Takahashi 
et al. (2009), the cool phase and hot phase coexist and dominate in the central 70/iy/ kpc 
of the Centaurus cluster, exhibiting a temperature ratio T c /Th « 0.46. These properties are 
analogous to those of A1795. In addition, a similar central excess in the total gravitating 
mass profile, with a radius of 50/ifj 1 kpc, has also been found in the Centaurus cluster (Ikebe 
et al. 1999). Thus, the two clusters are very similar in the ICM properties, as well as in the 
total mass distribution. 

As shown in §3.1.4 and Figure 46, the cool (~ 2 keV) phase component of A1795 
occupies up to only 20% volume even in the innermost region, and in Figure 6 and 7, no 
apparent separation on the spatial distribution is seen between the cool phase and hot phase 
ICM in the cluster core. This indicates that the cool phase is substantially mixed into the 
surrounding hot phase, instead of forming any "cool-phase-only" core on scale of > 15/if/ 
kpc. Since such a 2P structure is seen predominately around a cD galaxy, with an enhanced 
metallicity (particularly in the cool phase), the phenomenon is very likely related to the cD 
galaxy. 
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4.2. Bubble Uplifting Model 

One possible formation mechanism of the 2P structure could be an AGN-driven gas 
transport. As indicated by numerical simulations, e.g., Churazov et al. (2001) and Guo & 
Mathews (2010), the buoyant bubbles created by AGN outbursts can drag a certain amount 
of surrounding cool, metal-rich gas to larger radii. For example, the ICM in central 30/if/ 
kpc regions of the Hydra A cluster and M87 are known to be in a 2P or multi-phase form 
(Ikebe et al. 1997; Molendi 2002). More recently, the distributions of cool phase ICM 
were found to coincide with the powerful radio lobes and X-ray cavities, hence, the bubble 
uplifting model has successfully been invoked for both Hydra A cluster and M87 (Nulsen 
et al. 2002; Simionescu et al. 2008). These properties may be useful to explain the origin 
of the co-existing 2P ICM. As for A1795, by adopting the obtained distributions of total 
gravitating mass and cool phase gas mass (M(R) and M Cigas , respectively; §3.3.2), the energy 
required to uplift the cool component from cD galaxy center to current position (average 
R ps 40/ifj 1 kpc) is estimated as GM Cigas M(R') / R' 2 dR' ~ 10 58 ergs, which is consistent 
with the amount of mechanical energy injected from its central AGN outburst (Rafferty et 
al. 2006). Nevertheless, the faint X-ray cavity detected in A1795 is located < 20/if/ kpc 
from the center (Birzan et al. 2004), apparently much smaller than the radius of the cool 
phase ICM as shown in the obtained 2-D temperature and the Q c /Qh maps (Fig. 6 & 7). 
Hence the appearance of the 2P ICM cannot be solely ascribed to bubble uplifting via recent 
AGN activity. Neither can the cool phase ICM be bubble remnant, since the uplifted gas 
would conductively mix and reach thermal equilibrium with the surrounding hot phase ICM 
quickly, i.e., within ~ 10 6 yrs (e.g., Simionescu et al. 2008), after the bubble is torn apart. 



4.3. cD Corona Model 

The most natural and effective way to sustain a stable 2P structure is to separate 
the two ICM phases by magnetic fields. In such a case, the two phases become thermally 
insulated with each other, since the gyroradius of a thermal electron in a 1 — 20 fiG field 
is smaller than its mean free path by about 10 — 11 orders of magnitude (Sarazin 1988). 
Actually, cluster central regions are often threaded by rather strong (~ 10 /iG; Ge & Owen 
1993) magnetic fields, which are possibly related to the cD galaxy. Employing a general 
topological classification of such cD-galaxy-related magnetic field lines into closed and open 
ones, Makishima et al. (2001) proposed that the closed loops (with their both ends anchored 
to the cD galaxy) are filled with the cool phase ICM, while the open-line regions, connecting 
to outer regions, are permeated by the hot phase. Thus, the cool phase component is 
expected to emerge as numerous filamentary substructures in the core region, while the hot 
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phase distributes throughout the vast cluster volume. In addition, the cool phase is naturally 
more metal-rich than the hot phase in this scenario, as the supernova yield of the cD galaxy 
may be largely confined within the loops (Takahashi et al. 2009), which are surrounded by 
intruding, less contaminated hot phase ICM. This "cD corona" model is expected to provide 
a natural account for the appearance of metal-rich cool phase ICM in the center of A1795. 

As the cool phase is thermally insulated from the hot phase, and the radiative transport 
between the two phases can be ignored (Sarazin 1988), a continuous, efficient heating source 
located in the central region is required to prevent it from collapsing due to radiative cooling 
(see Aschwanden et al. 2001 for a review). As first pointed out by Rosner, Tucker & Vaiana 
(1978; hereafter RTV) for solar corona, and applied successfully to the Centaurus cluster by 
Takahashi et al. (2009), a loop structure as illustrated in Figure 11 has a built-in feedback 
mechanism to maintain thermal stability for the cool phase plasma. In order to characterize 
the loop heating mechanism in a quantitative way, we employed an updated version of RTV 
model, an analytical solution to the hydrostatic model introduced in Aschwanden & Schrijver 
(2002; AS02 hereafter) for solar corona and applied it to the 2P ICM. This model considers 
a thin, arch-like loop delineating magnetic field lines, for which the pressure between the 
loop-interior and surrounding plasmas are in equilibrium. To constrain the temperature 
and density structures of the loop, AS02 assumed mass and momentum conservations along 
the loop, as well as a total energy balance involving radiative loss rate, heating rate, and 
thermal conductive flux. The heating rate function has been defined in form of E R = E 
exp(—h/hn), where E is heating rate at the bottom of loop, h is the height in the loop 
plane, and h H is the scale height of the heating source (Fig. 11). The analytic solution 
by AS02 gave a loop temperature distribution increasing with h in general, and the exact 
temperature gradient was predicted to be sensitive to the heating scale hn. As shown in 
Figure 9 of AS02, for a uniform heating, i.e., > H, where H is the loop height, the loop 
temperature exhibits a mild decrease towards the loop bottom by > 60%, while for a small 
scale heating at the bottom, i.e., hu «C H, a large part of the loop becomes near isothermal, 
with only a sharp, small-scale temperature drop seen at the footpoints. Such a loop also 
obey a scaling law, as defined in Eq.(29) of AS02 (also see Eq.19 below), which describes the 
loop maximum temperature T max as a function of hu, H, and external pressure at footpoint 
Pq. Another important property of the loop-confined plasma is that, as shown in Eq.(30) of 
AS02, T max shows rather weak dependence on the heating rate E , i.e., T max oc E^ 7 , because 
a decrease (or increase) in Eq will make the loop thinner (or thicker), so as to keep the 
plasma luminosity nearly equal to Eq. 

Specifically, following AS02, the temperature distribution of the loop- interior plasma is 
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is directly derived from the scaling law (Eq.29 of AS02), g is the gravitational acceleration 
calculated from the obtained total mass profile (§3.3.2), a(H/hn) and b(H /hn) are the best- 
fit analytic functions, as defined in Eqs.(20) and (21) of AS02, respectively, to reproduce the 
dependence of T(h) on H/hn derived in a numerical way, and similarly, S± s (H, fin) is the 
best-fit empirical function of scaling law factor given in Eq.(35) of AS02. We adopted the 
best-fit coefficients of the hydrostatic solutions presented in Table 1 of AS02 to characterize 
a(H/hii), b(H/hu), and Sf s (H, /i H ), and then, to calculate T max , T(h)/T max , and n(h) based 
on the loop properties. In short, given a set of (Pq, H, hu), we may determine the temperature 
and density distributions of the loop-interior plasma. 

Next we applied the AS02 model to the X-ray observations of A1795 and, for a compar- 
ison, the Centaurus cluster, both harboring a prominent cool phase in the central region. In 
the calculation we determined the external pressures P by the temperatures and densities of 
the hot phases measured at the innermost regions, and took the loop heights H ~ kpc 
and 70h^ kpc for A1795 and the Centaurus cluster, respectively, as indicated by the radii 
of cool phases (§3.1.2; Takahashi et al. 2009). Figure 12a shows the predictions of Eq.(19), 
i.e., T max as a function of fin, in comparison with the measured maximum cool phase tem- 
perature, i.e., T max 2.4 keV and 2.2 keV for A1795 and the Centaurus cluster (Fig. 4a; 
Takahashi et al. 2009), respectively. This comparison yields h H ps 11/if/ kpc for A1795 and 
18hji kpc for the Centaurus cluster, which indicates that the heating is concentrated at the 
bottoms of the loops. 

Since the sizes of radio lobes of the central AGNs in A 1795 and the Centaurus cluster 
(« lOhyi kpc and 15/if/ kpc, respectively; Ge & Owen 1993; Taylor et al. 2002) agree well 
with hn estimated above, we speculate that the AGN feedback is a candidate heating source 
for the coronal loops. The life time of the radio lobes can be estimated by t sync = 4 ^2 7 ~ 
10 7 — 10 8 yrs, roughly consistent with the cooling time of the cool phase ICM. Furthermore, 
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as shown in Rafferty et al. (2006), the heating rates provided by the central AGN outbursts 
in A1795 and the Centaurus cluster were estimated to be ~ 1.6 x 10 44 ergs s" 1 and 7.4 x 10 42 
ergs s -1 , respectively, which agree well with their cool phase luminosities (~ 1.4 x 10 44 ergs 
s _1 for A1795 and 1.0 x 10 43 ergs s _1 for the Centaurus cluster). 

Having obtained the estimates of H and hu for the two clusters, we then calculated, via 
Eq.(17), the T(/i)/T max profiles, and show the results in Figure 126. Thus, a mild inward 
temperature decrease is expected, especially for the Centaurus cluster. Actually, Takahashi 
et al. (2009) showed (in their Fig. 4a) that the cool phase temperature of the Centaurus 
cluster decreases mildly, from ~ 2.2 keV (as already used in Fig. 12a) at ~ 25/if/ kpc , 
to ~ 1.6 keV = 0.72 T max at < 7/if/ kpc. This latter value, when plotted on Figure 12b 
(with a dashed cross), agrees with the prediction. The same is true in the A1795 case, which 
shows a more flat T(h)/T max profile that is consistent with the measured value of ~ 0.94 at 
15/if/ kpc (solid cross on Fig. 126). As shown in Figure 12c, we applied the obtained T(h) 
to Eq.(18) and calculated the normalized loop gas density profile n(h)/no for A 1795, which 
again agrees well with the observed profile. 

In summary, the cD corona view, combined with the AS02 modeling of the loop-confined 
plasmas, can account for several important properties of inner regions of A1795 and the 
Centaurus cluster. These properties include the stable 2P structure, the higher metallicity 
of the cool phase, the absolute values of T c , and the spatial T c and density distributions. 

5. CONCLUSION 

By analyzing the Chandra, XMM-Newton, and Suzaku data of the X-ray bright galaxy 
cluster A1795, we report clear preference for the 2P ICM model in the central 80/iff 1 kpc, 
which consists of a cool phase (2.0 — 2.2 keV) and a hot phase (5.0 — 5.7 keV) component. 
This 2P model provides significantly better fit to the deprojected spectra than the IP model 
with continuous temperature profile, while the latter cannot be fully ruled out based on 
current data. Combining the Suzaku XIS and the XMM-Newton EPIC & RGS, we have 
marginally detected a third weak 0.8 keV component in the inner 144/iyj 1 kpc region that 
can be ascribed to a portion of ISM component of the cD galaxy. Based on a 2-D spectral 
analysis with the ACIS data, we have revealed a possible spatial correlation between the 
cool phase and metal-rich gas in the 50 — 100/if/ kpc region. A follow-up XIS analysis shows 
consistent result, that the cool phase ICM does exhibit a higher metallicity (~ 0.80 Z & ) than 
the hot phase one (~ 0.36 Z & ). Hence, we have successfully resolved a 2 keV, metal-rich 
component associated with the cD galaxy, which is spatially mixed but thermally separated 
with the surrounding 5 keV cluster component. All these properties can be explained by the 
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cD corona model incorporating with the AS02 solution for quiescent coronal loops. 
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Table 1. Summary of Chandra, XMM-Newton, and Suzaku Observations 



Date 


Detector 


UDSliJ 




Dec 


Raw/Clean Exposure 


dd mm yyyy 






(h m s; J2000) 


(d m s; J2000) 


(ks) 


20/12/1999 


Chandra ACIS-S 


494 


13 48 56.5 


+26 36 26.1 


19.8/17.5 


21/03/2000 


Chandra ACIS-S 


493 


13 48 49.2 


+26 36 27.3 


19.9/18.9 


10/06/2002 


Chandra ACIS-S 


3666 


13 48 48.9 


+26 34 32.3 


14.6/13.7 


14/01/2004 


Chandra ACIS-S 


5287 


13 48 55.0 


+26 36 35.0 


14.5/13.5 


18/01/2004 


Chandra ACIS-I 


5289 


13 48 55.1 


+26 36 45.0 


15.2/14.7 


23/01/2004 


Chandra ACIS-I 


5290 


13 49 00.8 


+26 42 07.5 


15.1/14.7 


20/03/2005 


Chandra ACIS-I 


6159 


13 48 32.9 


+26 40 45.4 


15.1/14.6 


20/03/2005 


Chandra ACIS-S 


6160 


13 48 49.4 


+26 36 27.1 


15.0/14.1 


28/03/2005 


Chandra ACIS-I 


6161 


13 49 19.5 


+26 31 05.4 


13.8/13.3 


28/03/2005 


Chandra ACIS-I 


6162 


13 48 48.0 


+26 36 21.7 


13.8/13.3 


31/03/2005 


Chandra ACIS-I 


6163 


13 48 47.6 


+26 36 16.4 


15.1/14.8 


26/06/2000 


XMM-Newton EPIC-pn 


0097820101 


13 48 52.3 


+26 35 21.6 


47.7/30.2 


26/06/2000 


XMM-Newton EPIC-MOS1 


0097820101 


13 48 52.3 


+26 35 21.6 


50.1/43.6 


26/06/2000 


XMM-Newton EPIC-MOS2 


0097820101 


13 48 52.3 


+26 35 21.6 


50.1/40.3 


26/06/2000 


XMM-Newton RGS 


0097820101 


13 48 52.3 


+26 35 21.6 


66.6/40.1 


13/01/2003 


XMM-Newton EPIC-pn 


0109070201 


13 48 40.0 


+26 22 16.4 


70.3/48.1 


13/01/2003 


XMM-Newton EPIC-MOS1 


0109070201 


13 48 40.0 


+26 22 16.4 


64.4/54.1 


13/01/2003 


XMM-Newton EPIC-MOS2 


0109070201 


13 48 40.0 


+26 22 16.4 


64.4/54.2 


10/12/2005 


Suzaku XIS 


800012010 


13 48 53.8 


+26 36 03.6 


13.1/13.0 


10/12/2005 


Suzaku XIS 


800012020 


13 48 53.3 


+26 47 57.5 


24.4/24.0 


11/12/2005 


Suzaku XIS 


800012030 


13 48 53.5 


+26 59 58.2 


30.6/ a 


11/12/2005 


Suzaku XIS 


800012040 


13 48 53.5 


+26 24 02.5 


26.1/25.5 


12/12/2005 


Suzaku XIS 


800012050 


13 48 53.5 


+26 12 00.4 


40.1/37.2 


10/12/2005 


Suzafcu HXD 


800012010 


13 48 53.8 


+26 36 03.6 


10.4/10.4 



a This dataset is not used in analysis. See text §2.1.3. 



Table 2. Best-fit Temperature and Metal Abundance Gradients Obtained with 1P/2P 

APEC Fittings a 



Region /i C m 


Tip 


A 1P 




T c 


T h 


A 2 p 




(/if/ kpc) 


(keV) 


(Z Q ) 




(keV) 


(keV) 


(Zq) 





Chandra ACIS 






- 30 


99.8% 


q 97 +U.Z/ 

Ci ' z ' -0.31 


o.64±r,:oi! 


30 


- 51 


99.6% 


o 01 +0.27 
°- ol -0.22 


o.84l°; ° 


51 


- 80 


99.0% 


4 fi7 +0-26 
4 - ' -0.27 


0.6018;" 


80 - 


- 116 


97.8% 


5 61+ ' 37 
°- Di -0.38 


56+ 08 


116 


- 238 


93.4% 


6 "S2+ ' 32 


n oc+0.11 

0.36I 10 


238 


- 354 


86.7% 


a -in+0.34 
6 - 19 -0.39 


u "-' o -0.10 


354 


- 707 


82.4% 


r q-,+0.52 


0.281°;°? 


707- 


- 1335 


10.2% 


q 42 +0.77 
°-^ z -0.62 


0.23l°;i 2 2 



1750/1496 



2.18 
2.51 



+U.4b 



-0.18 
,+0.55 
-0.31 
9 19 +0.40 
z - ±z -0.44 



6.28 
5 

5.99 



+ 1.86 

1.80 
+0.92 
-1.00 
+0.64 

0.84 



0.54 
0.74 
0.51 



+0.11 
0.12 
+0.13 
0.12 
+ 0.09 
0.11 



1608/1494 



XMM-Newton EPIC 



0- 


- 30 


99.8% 


q n7 +u.ov 


52+ U '" b 
u - oz -0.05 


1040/881 


2 fi4 +u ' a '" 1 
z - u ^-0.64 


4.78t?& 

4 7<-|+0.65 
1- '>J-0.61 


n 45+0-06 
u - 4o -0.06 


929/879 


30 


- 51 


99.7% 


q f.y+0.06 

^•"'-o.oe 


61 +°' 06 
u - Di -0.06 




2 00+ ' 41 
z - uu -0.55 


n tro + 0.06 
U - b °-0.05 




51 


- 80 


99.5% 


4 16+ 008 
4 -!O_0.08 


o 5n+ ' 07 




i oq+0.53 
l.oo_ 3g 


5 in+o- 56 

J - lu -0.3O 


o 5i+o- 05 

u - oi -0.05 




80 - 


- 116 


99.1% 


5 01+ ' 12 


o 5n +0 ' 04 












116 


- 238 


97.3% 


t- q 9 + 0.10 

u -°^-o.io 


u - 4u -0.07 












238 


- 354 


91.7% 


c 7S +0.19 
a -'°-0.19 


n r, 7 +0.06 
u - z '-0.05 












354 


- 707 


85.5% 


t; O1+0.32 
°- zi -0.36 


25+ ' 05 
u - zo -0.06 












707- 


- 1335 


13.8% 


2 Q2+0.45 
z - yz -0.52 


n 20+°- 10 

u - zu -0.09 













Suzaku XIS 



- 


144 


99.4% 


4.69 


- 


144 b 


99.4% 




144 - 


- 320 


96.4% 


5.29 


320 - 


- 700 


86.9% 


5.39 


700 - 


- 1444 


14.3% 


3.20 



+ U.(lb 

-0.06 



+0.10 
-0.10 
+0.38 
-0.37 
+0.44 
-0.43 



0.48 



+ D.II3 

-0.03 



845/775 



0.36 
0.28 
0.23 



+0.04 
0.03 
+0.05 
0.04 
+0.12 
0.11 



2.09 
2.34 
[2.09 



+ U.4H 

0.45 
0.55 
0.84 
+0.49, 
0.451 



5 
5.58 
[5.48 



+0.41 

-0.36 
+0.51 
0.46 
+0.41, 
0.361 



0.47 
0.48 
0.35 



+0.03 
0.03 
0.03 
0.03 
+0.03 
0.03 



802/773 
791/771 



a Fraction of ICM emission in the total 0.5 — 10.0 keV spectrum(col. [2]), and the best-fit gas temperatures and 
metal abundance gradients obtained by applying IP ICM model (col. [3] & [4]) and 2P ICM model (col. [6] — [8]) to 
the ACIS, EPIC, and XIS dcprojcctcd spectra. The IP ICM and 2P ICM models are defined in §3.1.1 and §3.1.2, 
respectively. \ l v of the fits are shown in col. [5] & [9]. 

b Best-fit 2P gas temperatures and abundances, obtained by applying 2P ICM plus 0.8 keV component model 
(§3.1.3) to the XIS spectra extracted in — 144/i^j 1 kpc region. 

c The 2P temperatures in this region are fixed to the best-fit values in — 144/i^j 1 kpc region. 
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Table 3. Best-fit Metal Abundance and Absorption Gradients Obtained with the 2P 

VAPEC Fittings a 



Region 


O 


Mg 


Si 


Fe 


^H.Suzaku 


X 2 /" 


Nh, Chandra 


^H,XMM 


(hyi kpc) 


(2©) 


(Ze) 


(Z©) 


(Z©) 


(10 20 cm- 2 ) 




(10 20 cm- 2 ) 


(10 20 cm- 2 ) 


- 144 
144 - 320 
320 - 700 


u - 44 -0.22 

0.46^ 
0.54+°; 42 


u.oo_ 24 

50+ - 27 
u.ou_ 26 

n qq+o.52 


0-37±g;ig 
0.l0tg;» 


n 40+O.O3 
U - 48 -0.03 

n oi+0.05 
"••"-O.OB 


3 26 +0-54 

2 n7+°- 81 

1 41+0.51 


343/328 
299/284 
145/149 


9 , 7 +0.Bl 
-0.49 

1 20+ - 24 
1 - zu -0.32 

1 21+0- 25 


, no + 0.33 
1 -34+0.26 

1 n7+°- 25 

i - u ' -0.25 



a Best-fit metal abundance (col. [2]— col. [5]) and absorption gradients (col. [6]), obtained by applying the 2P ICM model 
to the deprojected Suzaku XIS spectra extracted in — 144h^ kpc and 144 — 320h^ kpc, and by applying the IP model to 
the XIS spectra in 320 — 700/iJ-j 1 kpc (see text §3.1.6). x 2 ' jv of the fits are shown in col. [7]. As a comparison, in col. [8] & 
[9] we show best-fit absorptions obtained with the deprojected Chandra ACIS and XMM-Newton EPIC spectra, respectively. 





Table 4. 


Analysis of X-ray Surface Brightness Profiles 






Projected and Double-/3 Fits a 




fi Model 






Double-/? Model 






Energy (kcV) 


Re" 




flcl 01 


Rc2 


02 


X z l» 


0.5 - 3.0 
3.0 — 8.0 


48.4 ± 0.9 0.56 ±0.01 
60.6 ± 1.2 0.60 ± 0.01 


810/143 
484/121 


53.2 ±2.1 0.72 ± 0.03 
53.5 ± 2.6 0.76 ± 0.04 


203.0 ± 7.2 
185.5 ± 16.5 


0.73 ± 0.05 
0.81 ± 0.08 


143/140 
135/118 






Deprojected IP (i, IP NFW, and 2P /3 + Double-0 Fits' 












IP Model d 




IP NFW Modcl d 






Energy (kcV) 


Rc 







c 






0.5 - 3.0 
3.0 - 8.0 


110.0 ± 0.4 
108.0 ± 0.5 


0.67 ± 0.01 
0.70 ± 0.02 


734/143 281.0 ±0.3 
364/121 289.0 ± 0.3 
2P /3 + Double-/3 Model 


5.12 ± 0.10 
5.12 ± 0.12 


351/143 
286/121 




Energy (kcV) 




0c 


Rc.hl 0hl R c ,h2 


0h2 






0.5 - 3.0 
3.0 - 8.0 
0.5 - 8.0 


21.2 ± 4.64 

[25.1]° 
25.1 ± 3.32 


0.60 ± 0.10 

[0.58] e 
0.58 ± 0.09 


73.1 ±2.2 0.81 ± 0.05 274.0 ± 8.6 
68.3 ±3.1 0.84 ± 0.06 211.5 ± 8.4 
72.3 ±1.9 0.81 ± 0.05 275.5 ± 8.7 


0.81 ± 0.05 
0.84 ± 0.06 
0.83 ± 0.06 


146/137 
143/115 
161/151 





a Models and corresponding parameters are defined in §3.3.1. 

b Best-fit R c /Rc,l/Rc,2/Rs/Rc,c/R c ,hi/Rc,h2 ar e given in the unit of kpc. 

c Modcls and corresponding parameters are defined in §3.3.1. 

d Surface brightness profiles were extracted from > lOO/i - ^ 1 kpc regions. 

e Due to the insufficient cool phase counts in 3.0 — 8.0 keV, the cool phase /3 parameters were fixed to the best-fit 0.5 — 8.0 
keV results. 
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R.A. (J2000) 




Fig. 1. — Exposure-corrected 0.7 — 8.0 keV images of A1795 obtained with (a) the Chandra 
ACIS and (b) the Suzaku XIS, where the thin and thick annuli used in the respective spectral 
analysis are also plotted. 1' corresponds to about 72hj 1 kpc. 
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Fig. 2.— (a) NXB-subtracted Chandra ACIS-I (upper), XMM-Newton EPIC-pn (middle), 
and Suzaku XIS-FI (lower) spectra extracted from > 1200/if/ kpc regions that were used to 
create the background template (§2.2), plotted against the best-fit models composed of CXB 
(blue), Galactic (green), and ICM (red) components, (b) NXB-subtracted deprojected Chan- 
dra ACIS-I (upper), XMM-Newton EPIC-pn (middle), and Suzaku XIS-FI (lower) spectra 
extracted from 320 — 700/iy] 1 kpc, along with the corresponding best-fit CXB (blue), Galactic 
(green), and ICM (red) models, (c) NXB-subtracted deprojected Chandra ACIS-S (upper), 
XMM-Newton EPIC-pn (middle), and Suzaku XIS-FI (lower) spectra extracted from the 
central 80/ifj 1 kpc, kpc, and kpc regions, respectively, fitted with IP (red 

solid; §3.1.1) and 2P (blue solid; §3.1.2) models. The prediction of each best-fit 2P compo- 
nent is shown with a blue dash line. For clarity the CXB and Galactic emission models are 
omitted from the plots. 
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Fig. 3. — Deprojected profiles of (a) gas temperature and (b) abundance, obtained with the 
Chandra ACIS (red), XMM-Newton EPIC (green), and Suzaku XIS (black) data under the 
IP ICM assumption (§3.1.1). The abundances of the central bins were also studied with the 
2P ICM model (§3.1.2), and the results are plotted in diamonds. 



-39 - 




Fig. 4. — (a) Deprojected 2P temperature profiles obtained with the Chandra ACIS (red), 
XMM-Newton EPIC (green), and Suzaku XIS (black) data (§3.1.2). For comparison, the IP 
result is also plotted with the solid line, which was calculated as the best-fit to the ACIS 
IP temperature profile (Fig. 3a) using the analytic formula in Vikhlinin et al. (2006; Eq.6 
therein), (b) Volume filling factor of the cool component, rj c , derived by the deprojected 
analysis of the Chandra ACIS (red), XMM-Newton EPIC (green), and Suzaku XIS (black) 
spectra (§3.1.4). The result obtained by fitting the ACIS surface brightness in 0.5 — 8.0 keV 
(§3.3.1) is plotted with a solid line. The ACIS and EPIC points on both panels are slightly 
displaced for clarity. 
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Fig. 5. — Emission measure distribution of different thermal gas components in the central 
144/17/ k P c 

region of A1795, obtained by fitting a multi-temperature model to the Suzaku 
XIS spectra after the projection effect was corrected (§3.1.5). Error bars were measured at 
the 68% confidence level. Dotted line shows the prediction of the coronal loop model. As 
described in §4.3, this model cannot be used to constrain the emission measures of the 0.8 
keV component or the hot ICM component. 
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R.A. (J2000) 

Fig. 6. — (a) A projected gas temperature map, (b) its 68% uncertainties, (c) a gas abundance 
map, and (d) its 68% uncertainties, all obtained with the IP spectral analysis of the ACIS 
data without incorporating deprojection (§3.2.1). In panels (a) and (c), superposed are the 
X-ray intensity contours calculated from the 0.5 — 8.0 keV data. All maps are centered on the 
X-ray peak of A1795. The cross marks the position of the X-ray cavity reported in Birzan 
et al. (2004). 
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Fig. 7. — A detailed comparison between (a) projected 2P abundance map and (b) Q c /Qh 
map of the central region, obtained with the 2P spectral analysis of the ACIS data without 
incorporating deprojection (§3.2.1). In panels (a) and (b), the X-ray intensity contours and 
the 2P abundance contours are superposed, respectively. The map center and the meaning 
of the cross is the same as in Figure 6. (c) A scatter plot between the values in panel (b) 
and those in panel (a) in 50 — kpc. Error bars are given at the 68% confidence level. 

For clarity, the abundances were averaged over 10% intervals in Q c /Qh- Solid line shows an 
analytic fit to the abundance-Q c /Qh relation, (d) Two point correlation function between 
the obtained 2P abundance map and the Q c /Qh map (black). A reference profile (red) is 
given after randomizing the two maps (see §3.2.3 for more details). 
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(f) 



Fig. 8. — Radial surface brightness profiles extracted in (a) 0.5 — 3.0 keV and (b) 3.0 — 
8.0 keV, obtained with the Chandra ACIS (0 - 1000/if/ kpc) and Suzaku XIS (1200 - 
1800/if/ kpc) after the instrumental normalizations are readjusted (§3.3.1). The best-fit 
/3+double— (3 models, based on the best-fit deprojected 2P results, are plotted with solid 
lines. For comparison, the best-fit /3 model and NFW model, based on IP results that can 
describe the surface brightness profiles extracted from > lOO/i^ 1 kpc regions, are also plotted 
with dash and dotted lines, respectively. 
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Fig. 9. — Total gravitating mass distributions of A1795, taking into account all directions 
(black solid line) , and northern half only (green solid line) , which were calculated based on 
the deprojected 2P temperature distribution and /3+double-/3 fitting of the gas density profile 
(§3.1.2 and §3.3.1, respectively). For comparison, the corresponding mass profiles obtained 
with the best-fit IP spectral analysis and double-/? gas density models (§3.1.1 and §3.3.2, 
respectively) are also plotted with red solid lines. Error bars are at the 68% confidence level. 
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Fig. 10. — NXB-subtracted Suzaku XIS-FI spectrum extracted from the central 320/if 1 1 kpc 
(black) and HXD-PIN (red) spectrum, compared with the CXB (green) and ICM (blue) 
components. The 2P ICM components were determined in the < 10 keV range, and ex- 
trapolated into the HXD-PIN range (see §3.4). The systematic NXB uncertainty of 2.0% is 
plotted in cyan. 
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Fig. 11. — A sketch of one of coronal loops that possibly exist in the envelope of the cD 
galaxy of A1795 (§4.3). H is the loop height and hu is the scale height of the heating source 
as shown in gray. 
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Fig. 12. — (a) Predicted loop maximum temperatures T max as a function of heating scale 
heights tin, calculated for A1795 (solid) and the Centaurus cluster (dashed). The horizontal 
lines mark the observed T max for these clusters, (b) Predicted T(/i)/T max profiles of the loops 
for A1795 (solid) and the Centaurus cluster (dashed), where h is the height in the loop plane. 
The crosses show the observed T(/i)/T max for the innermost shells presented in §3.1.2 and 
Takahashi et al. (2009). (c) Predicted n(h)/no profile of the loops for A1795 (solid). For 
comparison, data points with vertical error bars show the observed n(h)/no profile obtained 
with the 2P /3+double-/3 fitting of the ACIS surface brightness distribution (§3.3.1). 



